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We perform a detailed analysis of all existing data (CELLO, CLEO, BaBar) on the pion-photon 
transition form factor by means of light-cone sum rules in which we include the NLO QCD radiative 
corrections and the twist-four contributions. The NNLO radiative correction together with the 
twist-six contribution are also taken into account in terms of theoretical uncertainties. Keeping 
only the first two Gegenbauer coefficients a 2 and 04, we show that the la error ellipse of all data 
up to 9 GeV 2 greatly overlaps with the set of pion distribution amplitudes obtained from nonlocal 
QCD sum rules — within the range of uncertainties due to twist-four. This remains valid also for 
the projection of the la error ellipsoid on the (02, 04) plane when including a&. We argue that it is 
not possible to accommodate the high-Q 2 tail of the BaBar data with the same accuracy, despite 
opposite claims by other authors, and conclude that the BaBar data still pose a challenge to QCD. 

PACS numbers: 12.38. Lg, 12.38.Bx, 13.40.Gp, 11. 10. Hi 



I. INTRODUCTION 

The pion-photon transition 7* + 7* — > tt° is an ex- 
ample par excellence of an exclusive process that can be 
adequately described within QCD by virtue of collinear 
factorization, provided both photon virtualities are suffi- 
ciently large (for a review, see Then, the transition 
form factor, defined by the correlator of two electromag- 
netic currents 



d 4 ze-^'V(P) \T{j M (z)j v (0)} |0) 
ie„u a mi4F^{Q 2 ,q 2 ) (1) 



with Q 2 = —q\ > 0, q 2 = —q\ > 0, can be recast in 
leading- twist approximation in the convolution form [HQ 

F^'*(Q 2 , q 2 ) = j'dxTiQ 2 , q 2 , ^ 2 ,x) 

x V M(x,£) + 0{5 2 /Q*) , (2) 

where the pion decay constant is f„ = 132 MeV, S 2 
is the twist-four coupling, and where we assumed that 
the photon momenta are sufficiently large Q 2 , q 2 3> m 2 . 
This way, the quark-gluon sub-processes, encoded in the 
hard-scattering amplitude T, can be systematically com- 
puted order-by-order within QCD perturbation theory. 
The binding (nonpcrturbative) effects are separated out 
and absorbed into a universal pion distribution amplitude 

(DA) of twist-two (p {2 \x,p 2 ), defined first in Q. The 
variation of (x, /i 2 ) with the factorization scale is 
controlled by the Efrcmov-Radyushkin-Brodsky-Lepage 
(ERBL) renormalization-group evolution equation 
This implies that the pion DA can be expressed in terms 
of the Gegenbauer polynomials Cn^ 2 (2.T — 1) to read 



where ip n (x) = 6x{ 1 - x) Cl 1 2 {2x - 1 ) and s I 
6x(l — x) is the asymptotic pion DA @, Q. 

The coefficients a n (n ^ 2) have to be determined 
within some nonpcrturbative model or be extracted from 
the data, taking into account evolution effects to ac- 
count for their momentum-scale dependence. In most 
theoretical analyses they are derived from the moments 
(£, N )ir = Jq 1 dx(2x — 1) N ip£^ (x, /1 2 ) with the normaliza- 
tion condition J Q dxifn (x, /i 2 ) = 1. 

A process with two photons of high virtuality cannot 
be easily measured experimentally. Experimental infor- 
mation is mostly available for an asymmetric configura- 
tion in which one of the photons is quasi-on-mass-shcll 
with q 2 — > [B-0]- In that case, perturbative QCD can- 
not be directly applied at a consolidated level because 
the real photon is emitted at large distances and has, 
therefore, a hadronic content which demands the appli- 
cation of nonpcrturbative techniques. Such a method is 
provided by light-cone sum rules (LCSRs) Q (a recent 
application of this method to the pion-photon transition 
form factor is given in Q) which augments QCD per- 
turbation theory in approaching the real photon via a 
dispersion relation for F 7 7 w in the variable q 2 , keeping 
the large variable Q 2 fixed, viz., 



F 7 * 7 ^(g 2 , g 2 ) = / ds p< fj s 2 



S + T 



(4) 



where the physical spectral density p(Q 2 ,s) approaches 
at large s the perturbative one: 

p PT (Q 2 , s) = ilmF 7 * 7 ^ (Q 2 , -s - is) . (5) 



^ 2 Hx,fi 2 )=^ s (x) + 



7 ' 

*-^n=2A,... 



(a* 2 ) ipn(x) , (3) The resulting LCSR can be written as follows: 



2 



2^^, ( Q 2) = 01 / exp 



A/ 2 



(6) 



r 



with the spectral density p(Q ,x) = (Q 2 + s)p PT (Q 2 ,s) 
and the abbreviations s = xQ 2 /x and xo = Q 2 /(Q 2 + so). 
The first term in © stems from the hadronic content 
of a quasi- real photon at low s < so, while the second 
one resembles its point-like behavior at higher s > sq. 
The hadronic threshold in the vector-meson channel has 
the value s — 1.5 GeV 2 , whereas m p — 0.77 GcV [loj ). 
The Borcl parameter M 2 entering the LCSR in (jSJ) is not 
varied aiming to obtain the best stability of the LCSR. 
It is actually specified via M 2 — M 2 _ pt / (x)q2 from the 
two-point QCD sum rule for the p-meson, where the cor- 
responding Borel parameter is M|. pt G [0.5 -r- 0.8] GeV 2 . 
Here, (x)q2 is some average value of a; at a fixed scale 
Q 2 in the integration region for the first integral on the 
r.h.s. of Eq. ([5]), i.e., x (Q 2 ) < (x)q2 < 1 — see also [1]. 

The remainder of the paper is organized as follows. In 
the next section, we review the current understanding 
of the experimental data on the pion-photon transition 
form factor, focusing on the problems related to the un- 
expected rise of the BaBar data beyond 10 GeV 2 . In 
Sec. IIII1 we present the salient features of a calculation of 
the pion-photon transition form factor using LCSRs with 
NLO accuracy and including the twist-four contribution. 
The results of our calculation are presented and discussed 
in detail in Sec. lIVl in terms of two tables and several fig- 
ures, paying particular attention to the intrinsic theoreti- 
cal uncertainties stemming from various sources. Finally, 
a summary of our main findings is given in Sec. |Vj where 
we also draw our conclusions. 



II. STATUS OF EXPERIMENTAL DATA 

Using the method of LCSRs with fixed M 2 = 
0.7 GeV 2 , Schmedding and Yakovlev (SY) (TJ performed 
a next-to-leading order (NLO) QCD analysis of the data 
on F 7 * 77r , obtained by the CLEO Collaboration @, with 
the inclusion of twist-four contributions,which they mod- 
eled by the asymptotic form, notably, || [lj| , 



^\x^ 2 ) = {^/3)5 2 {^ 2 )x 2 {l-xf 



(7) 



0.56 and 0,4 = were found to be outside the 2a level. 1 
Later, this type of analysis was refined in [III, [l|| by 
taking into account the correct ERBL evolution of the 
pion DA to each measured momentum, including also 
the variation of the twist-four coupling S 2 = 8 2 {p\) = 
0.19 ± 0.04 GeV 2 in J7J). This way it was shown that the 
DAs v?| s and <pj z are outside the 3a and 4cr error ellipses, 
respectively, while the best fit yields x 2 df = 0-47 (where 
ndf denotes the number of degrees of freedom), see for 
more details in [r|. 

Even more important is the fact that the la error el- 
lipse 2 of the CLEO data strongly overlaps with the ad- 
missible values of the coefficients ai and CI4, determined 
before in the framework of QCD sum rules with non- 
local condensates (NLC SRs) [l7|. In particular, we 
proposed a "bunch" of admissible pion DAs with the 
central point af MS (/i 2 ) = 0.20 and a^ MS (pj) = -0.14 
(termed the BMS model) that turns out to be within 
the la error ellipse of the CLEO data [f|. These pion 
DAs (and similar ones obtained by two of us (BP) in 
[2 El by employing an extended version of the Gaus- 
sian model of the nonlocal QCD vacuum) have a dis- 
tinctive feature: their endpoints at x — and x = 1 
are strongly suppressed — even relative to the asymptotic 
pion DA. This suppression is due to the finite virtual- 
ity A 2 = {q{iga liv G^)q)/{2{qq)) = (0.4-0.5) GeV 2 of 

(2) 

the vacuum quarks, which enters the NLC SR for ip-x 
in terms of the various condensate contributions (more 
details can be found in [l7l. I20I]). 

These findings remain valid even if one assumes a much 
stronger contribution from the twist-four pion DA, as 
estimated in [l3| using the renormalon model 21 (see 
also [H). The calculation of F^*^ beyond the NLO 
was carried out in [23[ by taking into account that part 
of the next-to-next-to-leading order (NNLO) correction 
that is proportional to the first coefficient Bp of the 8 
function, known from the earlier work in [241 ] . It turns 
out that this radiative correction is negative, thus pro- 
viding suppression of the form factor at the level of about 
8%, slowly decreasing with increasing Q 2 . To the extent 
that one's ambition is to comply with the CLEO data on 
7*7 — > 7T°, it is mandatory to have a pion DA that has 



The major outcome of their analysis was that the 
best agreement with these high-precision data can be 
achieved by a pion DA with the coefficients af Y (^§Y = 
5.76 GeV 2 ) = 0.19, af Y (//§ Y ) = -0.14, which after evo- 
lution to the reference scale p 2 = 1 GeV 2 [l3| become 
of Y (p 2 ) = 0.27, af Y (/i 2 ) = -0.22. 

At the same time, the asymptotic pion DA and the 
Chernyak-Zhitnitsky (CZ) model [3 with a^ z (/i 2 ) = 



1 The value of aSp 2 was originally determined by Chernyak and 
Zhitnitsky at the (low) scale (j.q z = 0-25 GeV 2 Il4l and reads 
a^ z ((j,Q Z ) = 2/3. Its determination at the usual reference scale 
/if = 1 GeV 2 is described in App. B of Ref. [T5I ] ■ 

2 Here and below we denote by la ellipse (ellipsoid) a confidence 
region in the (02, £14) plane that represents a coverage probability 
equal to 68.27%. 
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its endpoints x = 0, 1 suppressed with a moderate dip 
at x = 1/2, just as the BMS model extracted from NLC 
SRs. It is worth mentioning that this model gives rise to 
a scaled form factor Q 2 F~ I 77r (<3 2 ) that approaches the 
asymptotic limit a/2/tt, predicted by perturbative QCD, 
from below, while at low momenta it also yields reason- 
able agreement with the CELLO data Q. 3 Therefore, 
it came as a surprise that the data on 77* — > ir°, pre- 
sented by the BaBar Collaboration in 2009 0, indicate 
above approximately 10 GeV 2 a marked growth of the 
form factor with Q 2 (see upperpoints (labeled by open 
crosses) in Fig. [1] taken from |25|). In fact, the BaBar 
Collaboration fitted their data with the expression 



0.3 



(3/5)Q 2 F n (Q 2 ) (GeV) 



Q 2 F 1 77r 



10 GeV' 



(8) 



where /3 = 0.25 ±0.02, i.e., F 7 * 77r - Q~ 3/2 . This behav- 
ior of the form factor is incompatible with the collinear 
factorization and exceeds the asymptotic prediction con- 
siderably up to the highest measured momentum value 
just below 40 GeV 2 . Moreover, the high-Q 2 BaBar tail 
contradicts the CLEO data which are well-described by 
the dipole fit Q 



Q 2 F~ I ' 77r 



Q 2 + A 2 



(9) 



with A sa 776 McV. This scale behavior becomes more 
puzzling by comparing the 7r° BaBar data with the pre- 
dictions extracted by the BaBar Collaboration from their 
analysis of the two-photon 77 and 77' decays [25[ — Fig. 
[TJ While the direct 7r° data grow with Q 2 and necessi- 
tate a wide DA with strong contributions from the end- 
point regions x = 0, 1, the converted n — rf data are in 
perfect agreement with the prediction obtained with the 
endpoint-suppressed BMS pion DA (solid line in Fig. [TJ . 
To describe both sets of the BaBar data with the same 
accuracy looks like a catch-22-situation because it de- 
mands an impossible compromise for the structure of the 
pion DA (more on this and its implications later). 

Many theorists have attempted during the past two 
years to explain the rise of the BaBar data using various 
schemes and drawing incongrucnt conclusions. We may 
group these theoretical proposals into four categories: 

(i) Abandon collinear factorization 
An increase of the form factor with Q 2 can be achieved 
by employing a flat-top pion DA, as proposed by 
Radyushkin [2(| and by Polyakov (27j — which used, how- 
ever, different schemes. The key element of these ap- 
proaches is a nonvanishing pion DA at the endpoints 
x = 0, 1. Such a pion DA ascribes the rise of the large- 
Q 2 tail of the BaBar data to a logarithmic increase of 
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FIG. 1: The 77* — > |n) transition form factor multiplied by 

3/5 in comparison with Q 2 F 77 *' r (Q 2 ) 0. The dashed line 
indicates the asymptotic limit for Q 2 F" n ^(Q 2 ). The dotted, 
dashed-dotted, and solid curves show the predictions from [l5l . 

for the asymptotic, the CZ, and the BMS DA, respectively 
(figure taken from [25l]h 

the form factor induced by the flat-top DA. Note, inci- 
dentally, that as shown in |28| , a In Q 2 behavior can actu- 
ally be imitated by including just the first three Gegen- 
bauer polynomials without demanding the flatness of the 
pion DA. However, one has to pay a high price for this 
agreement: already the NLO contribution at the default 
normalization scale \j? = Q 2 turns out to be huge for the 
flat-top DA, forcing a much lower normalization scale be- 
low Aqcd and making it virtually impossible to include 
evolution effects which arc at the heart of QCD. 

(ii) Deny rise of the form factor 
Very recently, Agaev et al. Q argued that the trend of 
the high-Q 2 tail of the 7r° BaBar data can be reproduced 
within QCD using LCSRs. Their analysis of the form 
factor includes radiative corrections at the NLO level, 
twist-four contributions, and, for the first time, also a 
twist-six term which, though it has a positive sign, turns 
out to be small for their choice of the Borel parameter 
M 2 = 1.5 GeV 2 to provide any serious enhancement. 4 
The distinctive feature of their models for the pion DA is 
a large and positive coefficient 04 to which they attribute 
the rise of the form factor at high Q 2 , while still higher 
coefficients a n up to n = 12 have smaller values. The best 
agreement with the ir° BaBar data 0] is provided by a 
modified flat-like DA with 02 = 0.13, which has a large 
tail of Gegenbauer coefficients, while the other two con- 
sidered models are given by truncated expansions with 



3 Note that at low momenta close to Q 2 £s 1 GeV 2 , the twist- 
two contribution becomes comparable in size with the twist-four 4 This fact is emphasized by these authors as being rather crucial 
term. in selecting that large value of the Borel parameter. 
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a n up to n — 8 and provide form-factor predictions that 
almost scale with Q 2 , but having a normalization exceed- 
ing considerably the asymptotic limit of V^f^ . Even the 
model with a„ coefficients up to n = 12, which become 
successively operative with increasing Q 2 , cannot simu- 
late the increase of the form factor sufficiently well in 
comparison with the BaBar fit © used in Q — provided 
one takes this fit seriously. At the same time, as one may 
appreciate from Fig. [1] all these models overestimate the 
BaBar data [25[ on the 77* — > 77(77') transition form fac- 
tor significantly. Taking into account the NNLO radia- 
tive corrections, would somewhat reduce the magnitude 
of the form factor by a few % [23], but incorporating 
this effect would also shift the prediction further away 
from the 7*7 —> ir° BaBar data — the main goal. Thus, 
the conclusion drawn in [9| that one may achieve simul- 
taneous agreement with all data on F 1 77r using LCSRs 
within QCD seems rather perfunctory. At the same time, 
the fact that the twist-six term turns out to have a posi- 
tive sign renders Chernyak's recent claim in [2!| that one 
can describe all data using the CZ DA supplemented by 
a power correction 1 /Q 6 with a huge negative coefficient 
-(1.2 GeV 2 ) 2 rather speculative. 5 

(iii) Contextual explanations 

Various authors (e.g., [3(|) have presented calculations 
of Q 2 F 1 77r within the context of particular models — 
appealing, for instance, to the triangle anomaly — and 
have indeed obtained a logarithmic increase, emulating 
this way the large-Q 2 tail of the BaBar data. The prob- 
lem with such approaches is that they involve (nonper- 
turbative) context-dependent mass scales that cannot be 
derived from QCD and are specific for this particular 
process. 

(iv) fej^ factorization 

Some authors try to explain the increase of the form 
factor by retaining the partonic fej^ degrees of freedom 
unintegrated and employin g a modified fcj^ factorization 
approach, see for instance [3lH34| . However, it is well- 
known from previous research [35l . l36j that Sudakov ef- 
fects and fcj_-generated power corrections give rather sup- 
pression than enhancement, albeit self- reinforcement is 
also possible |37[. This subject is, however, outside the 
scope of the present investigation; it will be discussed in 
detail elsewhere. 



III. DATA ANALYSIS 

The strategy of our investigation is to use LCSRs in 
order to process all available experimental data on the 
pion-photon transition form factor from different exper- 
iments [10] j making some technical improvements rel- 



ative to our previous works in [13J, |l5|, [l6|, l23| , and ex- 
tending the analysis to the level of three Gegenbauer co- 
efficients 02,04,06- The improvements and the applied 
algorithms are explained in detail below. 

To compute the form factor, we employ the LCSR 
(HU). In order to resolve the p and uj meson reso- 
nances, we employ not a simple <5-function ansatz, but a 
finite- width Breit-Wigner formjl, [23|. Adopting Khod- 
jamirian's general arguments [8J, we relate the Borel 
parameter M 2 in the LCSR © to the Borel param- 
eter M2.pt entering the two-point QCD sum rule for 
the p meson (38|. Hence, instead of using the fixed 
value M 2 = M 2 = 0.7 GeV 2 , as in our previous works 
[H El H H [2| , we write M 2 (Q 2 ) = M 2 / (x) Q 2 , where 
the average (x)q2 is calculated by means of the Borel 
exponential exp \Q 2 x/(Mq x)] in This prescription 
amounts to M 2 (Q 2 ) G [0.7-=- 0.9] GeV 2 . 

The key ingredients of our data-analysis procedure are 
the following: 

(i) The NLO gluon radiative corrections in the spectral 
density are taken into account usin g a corrected expres- 
sion with respect to Eq. (3.12) in |23| which does not 
contain the error pointed out in Note that this er- 
ror does not affect the results of our previous calcula- 
tions with LCSRs in [U, [H [H [H. We adopt the so- 
called default rcnormalization-scale setting and identify 
the factorization and the renormalization scale with the 
large photon virtuality Q 2 . This avoids the appearance 
of (large) logarithms of these scales. 

(ii) We take into account the twist-four contribution, al- 
lowing for a significant variation of the parameter <5 2 = 
0.19 GeV 2 in the range 0.15 GeV 2 to 0.23 GeV 2 , referring 
the reader for a detailed discussion to our previous analy- 
sis in [l6j . As already mentioned, the use of a nonasymp- 
totic form for ip^ would not affect the results signifi- 



cantly [13j,l22j. The parameter S 2 (p 2 ) is evolved with p 2 
according to the one-loop renormalization-group equa- 
tion, whereas other perturbative corrections to the twist- 
four contribution are ignored — expected to be small, 
(iii) The evolution of the coefficients a n is taken into ac- 
count at the NLO level of accuracy using the QCD scale 



parameters Aqq D = 370 McV and A^ D = 304 MeV, 
consistent with the NLO estimate a s (M§) = 0.118 [Tol | . 
(iv) The inclusion of the NNLOh radiative correction to 
the form factor, calculated in [24j and used in the spectral 
density in [23| . as well as the twist-six term, computed 
for the first time in Q, are taken into account implic- 
itly within the theoretical uncertainties. Their explicit 
evaluation is relegated to a future dedicated investiga- 
tion. This treatment of the two contributions in terms 
of uncertainties is justified by the fact that for the av- 
erage value of M 2 (Q 2 ) ~ 0.75 GeV 2 the net result ap- 
pears to be small, decreasing with Q 2 from +0.005 at 
Q 2 = 1 GeV 2 — where the twist-six term prevails — down 
to -0.003 at Q 2 = 40 GeV 2 — where the NNLO correc- 
tion becomes stronger. Note that in the calculation of the 
NNLO term only the convolution of the hard-scattering 



5 However, other power corrections of different origin with nega- 
tive sign cannot be excluded, but can hardly be accommodated 
within QCD. 
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FIG. 2: (color online). 3D graphics of la error ellipsoids in the space spanned by the Gegenbauer coefficients a n (n = 2, 4, 6) of 
all data [a-Q on the pion-photon transition form factor in the range [1 -j- 9] GeV 2 (left panel) and [1 40] GeV 2 (right panel) 
processed with the help of LCSRs (details in the text). The projection of the la ellipsoid on the plane (02, 0,4,) is represented, 
in both panels, by the larger ellipse in red color. The smaller ellipse denotes the cross section of the ellipsoid with the (02, 04) 
plane. The shaded (green) rectangle encloses the region of those 02, (24 pairs that are allowed by NLC SRs [T3|, whereas its 
middle point (Jt) marks the BMS pion DA. All results are shown at the scale after NLO evolution. 



amplitude Tp with the asymptotic DA <p as is taken into 
account (see for more details in [HI). It is worth men- 
tioning that this behavior is sensitive to the choice of the 
Borel scale adopted in our work. Would we use instead 
the value M 2 = 1.5 ±0.5 GeV 2 , advocated for in @, then 
the twist-six term would be much smaller and the net re- 
sult would be everywhere negative and almost constant: 
w -0.004. 



IV. DISCUSSION OF RESULTS 

To be able to make precise estimates of the influence 
of the high-Q 2 tail of the BaBar data on the form fac- 
tors, we use two different data sets: one which collects 
all available data from three experiments for the Q 2 
region between 1 and 9 GeV 2 (termed for convenience 
"CLEO regime" ) and a second set which includes as well 
the high-Q 2 BaBar data up to about 40 GeV 2 . These 
results are displayed in Fig. [5J left and right panel, re- 
spectively. In practice, this analysis procedure allows us 
to understand the role and relative strength of the coeffi- 
cient ae in determining the characteristics of the la error 
ellipsoid of each data set selection, shown in this figure 
by its principal axes and corresponding ellipses. 

The middle point (0.17,-0.14,0.12) of the la error 
ellipsoid in the left panel of Fig. [3J resulting from ana- 
lyzing the data in the CLEO regime, corresponds to a 
rather good x„df ~ 0-4. This ellipsoid is stretched along 
the axis and has a sizable intersection with the (02, 04) 
plane, shown in the figure as a smaller ellipse inside the 
green rectangle. In contrast, the inclusion of the high- 



Q 2 BaBar data (right panel) forces the ellipsoid away 
from that plane to a new middle point (0.18, —0.17, 0.31) 
with a worse value x„df ~ 1- Observe that the projec- 
tion of the ellipsoid on the (02,04) plane remains almost 
unchanged, while its middle point ▼ moves slightly away 
from its previous location, albeit it still resides inside the 
region of the negative values of 04 . The fact of the matter 
is that both error ellipses mentioned above cover a quite 
large region of those values of the coefficients 02 and 04 
which have been determined in [llj with the help of NLC 
SRs — slanted (green) rectangle in Fig. rjj In particular, 
the position of the BMS pion DA in the (02, 04) plane is 
just inside both ellipses for the CLEO regime with respect 
to all data (left panel) and still inside the 2D projection 
when one takes into account also the BaBar data up to 
40 GeV 2 (right panel). 

To further quantify these considerations, we perform a 
more detailed data analysis with the focus on the (02, 04) 
plane and show the results in Fig. [3J The left panel 
refers to the data selection for the Q 2 values in the CLEO 
regime, while the right panel shows the analogous find- 
ings for the second case when we implement the analysis 
by including all BaBar data up to the highest measured 
value of Q 2 . To this end, we calculate the la error ellipses 
by allowing the parameter S 2 to vary by 20% around the 
value 0.19 GeV 2 . All ellipses obtained this way are then 
merged together into a single distorted la error ellipse 
shown in Fig. [3] for different situations. These are the 
following: The largest ellipse — dashed line in red color 
with the middle point T — represents the result of com- 
bining the projections on the plane (02,04) of the 3D 
data analysis. The smaller ellipse (solid blue line) shows 
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FIG. 3: (color online). Distorted la error ellipses of the F Y >7 * (Q 2 ) data in the range [1-^9] GeV 2 (left panel) and [1 + 40] GeV 2 
(right panel) from various experiments [51471. l25j using different data-analysis procedures. These ellipses (see for explanations 
in the text) are the result of unifying different ellipses which correspond to the twist-four parameter varying in the range 
S 2 = 0.15 + 0.23 GeV 2 . The slanted rectangle in green color encloses the area of a 2 and 04 values determined by NLC SRs [171 ] . 
with the BMS pion DA being marked by K. The middle points of the ellipses are also labeled (+ and T) together with the 
asymptotic DA (♦), the CZ DA (■), and Model III from [y] (A). The vertical lines indicate the range of an values determined 
by two lattice simulations: [3t| — dashed lines; [4^] — dashed-dotted (blue) lines. The Borel parameter M 2 (Q 2 ) is varied with 
Q 2 in both panels as explained in the text. All results are shown at the scale /x| Y - 



the analogous result of a 2D analysis by means of 02 and 
04. Its middle point + with the coordinates (0.15, —0.09) 
and Xndi ~ 0.5, almost coincides with the middle point 
* of the area of values determined by NLC SRs. Finally, 
the smallest ellipse (thick line), entirely enclosed by the 
previous one, is obtained by combining the intersections 
with the (02,04) plane of all ellipsoids generated by the 
variation around the central value of 5 2 . The ellipsoid 
which corresponds to the central value S 2 = 0.19 GeV 2 
is shown in the left panel of Fig. [21 

The locations in the (02,04) plane of some character- 
istic pion DAs are also indicated in Fig. [3] notably, the 
asymptotic DA (♦), the CZ model (■), and also the pro- 
jection of Model III from Q (a). It is worth empha- 
sizing that the slanted (green) rectangle that contains 
those values of 02 and 04 derived from NLC SRs 17 1 



lies almost completely within both larger error ellipses — 
greatly overlapping with the smallest one as well. In 
particular, the BMS model DA X turns out to be inside 
of all la error ellipses. Hence, the results of the 2D and 
3D data analyses are in good mutual agreement, while 
we observe no tension between them and the theoretical 
predictions obtained from NLC SRs [It} — at least at the 
level of Xndf < 0.5. What is more, all calculated la er- 
ror ellipses are complying fairly well with the boundaries 
for 02 provided by two independent lattice simulations. 
The older estimate from [39| is indicated in both panels 
of this figure by the vertical dashed lines, whereas the 
very recent results given in Ref. (Table XII there) 
are denoted by the dashed-dotted vertical lines in blue 
color. 

As one can see from the right panel of Fig. [31 the situ- 
ation changes considerably when the high-Q 2 tail of the 
BaBar data [3] is included in the data analysis. Using the 
same designations as in the left panel, we show the anal- 
ogous unified error ellipses, observing that now the ellip- 



soid has no intersection with the (02,04) plane. More- 
over, the composed error ellipse resulting from the 2D 
analysis (solid blue line) moves out of the region of the 
negative values of 04 and inside its positive domain. This 
is accompanied by the significantly worse value Xndf ~ 2, 
in sharp contrast to the value Xndf ~ 0-5 we found for 
the data set in the CLEO-regime. On the other hand, 
the unified la error ellipse of the 3D projections on the 
(02,04) plane (larger dashed red ellipse), keeps its po- 
sition unchanged and still harbors a big portion of the 
area for the 02, 04 values, compatible with the results 
from NLC SRs (shaded rectangle in green color). 

Comparison with the lattice findings shows that the 
3D error ellipse lies almost entirely within the boundaries 
from [3^| (dashed vertical lines) but intersects with the 
interval determined in (dashed-dotted vertical lines) 
only for the larger values of 02. The analogous ellipse of 
the 2D analysis only poorly complies with this small 02 
window of [40j, whereas it partly overlaps with the low 
end of the 02 range determined in [39j . Thus, the mutual 
agreement between the 2D and the 3D analysis, observed 
in the left panel, now deteriorates. Moreover, it becomes 
obvious from this figure that Model III from has a 
projection on the (02,04) plane that lies outside of all 
considered la error ellipses of the data. Note incidentally 
that selecting for the Borel parameter the value M 2 = 
1.5 GeV 2 , as used in Q, we find for their Model III an 
agreement with the BaBar data on the level of x 2 df > 1.5. 

This discussion becomes more substantiated by in- 
specting Fig. [¥] which compares theoretical predictions 
with various experimental data with appropriate desig- 
nations provided in the figure. Our calculations are dis- 
played in the form of a shaded strip (green color) with a 
width reflecting various theoretical uncertainties. These 
are (i) uncertainties owin g to the spread of the "bunch" of 
the pion DAs derived in [17] . (ii) those originating from 
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TABLE I: Results for the scaled pion-photon transition form factor calculated in this work (last column) in comparison with the 
values measured at the same momentum scales by the CLEO 0| and the BaBar [3, EEf Collaborations, both for the 7*7 — > rr 
and the 7*7 — ¥ |n) decays. The last entry at the timelike momentum Q 2 = 112 GeV 2 is taken from [4l[. The origin of the 
theoretical uncertainties in the last column is explained in the text. 



CLEO/BaBar 


CLEO/BaBar 


CLEO 


BaBar 


BaBar (77, 77') 


This work 


Q 2 interval 


Q 


Q I^CLEO (Q )l 


Q I^BaBar (Q ) 1 


(3/5)Q 2 |F^ ar (Q 2 )| 




[GeV 2 ] 


[GeV 2 ] 


[0.01 x GeV] 


[0.01 x GeV] 


[0.01 x GeV] 


[0.01 x GeV] 


1.5 


- 1.8 


1.64 


12.1 ±0.8 ±0.3 


— 


— 


11 1 A'^'W T U -° 

11 ' 1 ^-1.02-0.3 


1.8 


- 2.0 


1.90 


11.7±0.7±0.3 


— 


— 


„ . _i_ i nn 1 ft R 

1 ri/i -r-L-uy +U.O 
iZ ' U4 -1.03-0.3 


2.0 


- 2.2 


2.10 


13.8 ±0.8 ±0.3 


— 


— 


1 1 hq _i_n ^ 
10 CQ~r 1 - l Jy -|-U.O 

±z - oo -1.01 -0.3 


2.2 


- 2.4 


2.30 


12.7 ±0.9 ±0.3 


— 


— 


1 n^i^ -f-u.o 

iJ ' uo -0 .99-0 .3 


2.4 


- 2.6 


2.50 


13.5 ±1.0 ±0.3 






H „ . _|_1 C[A _L_0 ^ 

10 '^ -0.95-0.3 


2.6 


- 2.8 


2.70 


15.1 ± 1.1 ±0.4 






., rt r-,t -1-1 01 -1-0 5 

lo. 1 r_ .92 -0.3 


2.8 


- 3.1 


2.94 


13.7 ±1.2 ±0.3 


— 


— 


. nn _i_n no 1 n r 
i4,uz -0.90-0.3 


3.1 


- 3.5 


3.29 


14.5 ±1.2 ±0.4 


— 


— 


1 a ne 1 n c 

14.38±°;l|±°;i 


3.5 


- 4.0 


3.74 


13.2 ±1.4 ±0.3 


— 


— 


, i_a no _i_n K 

^■'^-0.85-0.3 


4.0 


- 4.5 


4.24 


13.4 ±1.5 ±0.3 


15.04 ± 0.39 


— 


±O.U^_ .S2-0.3 


4.0 


- 5.0 


4.44 






14.89 ±0.26 


±,J - lz _0.80-0.3 


4.5 


- 5.0 


4.74 


15.4 ±1.7 ±0.4 


14.91 ±0.41 




■ _ „__|_o S3 -1-0 5 
±O.ZO_q 79 _ 3 


5.0 


- 5.5 


5.24 


14.5 ±1.8 ±0.4 


15.74 ±0.39 


— 


. _ . n on _i_n ^ 
LO.W_ 76 _Q 3 


5.0 


- 6.0 


5.45 


— 


— 


14.86 ±0.27 


1 a 7n 1 n c 
lO.^»_ 75 _ 3 


5.5 


- 6.0 


5.74 


15.5 ±2.2 ±0.4 


15.60 ±0.45 


— 


1 ^ ^s^ u- ' ' + u - 3 

±O.OO_ .74_o.3 


6.0 


- 7.0 


6.47 


14.8 ±2.0 ±0.4 


16.35 ±0.36 


— 


, .in 7A±fl R 

-L-J- ' "_o.71 -0.3 


6.0 


- 8.0 


6.85 






15.52 ±0.28 


_ _ „ . _i_a 70 in e 

ID. o^_ ,7o_o.3 


7.0 


- 8.0 


7.47 




16.06 ±0.47 




. _ ,-.,.-4-0 71 -1-0 5 


7.0 


- 9.0 


7.90 


16.7 ±2.5 ±0.4 


— 


— 


1 a en 1 n K 
■t AQ+U.09+U.0 

1U.UO_ .67-0.3 


8.0 


- 9.0 


8.48 


— 


16.73 ±0.60 


— 


, j-, ,-,4-0 fiS -1-0 R 

Ifi 1 QtU.OO tU.O 

±U.±^_ .65_o.3 


8.0 


- 10.0 


8.87 


— 


— 


16.27 ±0.37 


1 A *7 [AC 

1 A 1 7+0.67+0.5 
J-"- 1 ' -0.65 -0.3 


9.0 


- 10.0 


9.48 


— 


18.53 ±0.55 


— 


1 a (i k 1 n R 
1 A 0'\"'" u ' Da -f-u-o 
lO. ZO_ 64_o.3 


10.0 


- 11.0 


10.48 




18.66 ±0.76 




iU ' JU -0 .62-0.3 


10.0 


- 12.0 


10.90 




— 


16.76 ±0.47 


1 A Q 1 A c 

-1 f} A +O.DO +0.5 

J- -^ -0.61-0.3 


11.0 


- 12.0 


11.49 


— 


19.16 ±0.78 


— 


1 A SO 1 A K 

1B -*°-0 .61-0 .3 


12.0 


- 13.5 


12.71 




17.5 ±1.1 




lD.dO_ 59 _ 3 


12.0 


- 14.0 


12.91 






15 95 ± 65 


IK C7+0-60+0.5 
ru.O ( _q.59 —0.3 


13.5 


- 15.0 


14.22 




19.8 ±1.2 




16.65±°:£i°;« 


14.0 


- 17.0 


15.31 






17.83 ±0.67 


IK 79+0.57+0.5 
ID. '^-0.57-0.3 


15.0 


- 17.0 


15.95 




20.8 ±1.2 




IK 7c;+0.57+0.5 
ID. ' O— 0.56 —0.3 


17.0 


- 20.0 


18.40 




22.0 ±1.3 


17.06 ±0.97 


1 R oc+0.55+0.5 
1D.8D_q 55 _ 3 


20.0 


- 25.0 


22.28 




24.5 ±1.8 


17.46 ±1.10 


16 qq+o.53+0.5 
±u.»y_ 53 _ 3 


25.0 


- 30.0 


27.31 




18.111:2 


16.98 ±1.94 


17 1 1+0.51+0.5 
J -'- 1J --0.51-0.3 


30.0 


- 40.0 


34.36 




28.5±|;| 


17.95 ±2.26 


17 OO+0.49+0.5 
1 1 •^•J-0.49-0.3 






112 






20.16 ±1.65 


17 fio+0.42+0.5 
1 I .UO_ 43 _ 3 



the variation of the twist- four parameter S 2 , and (iii) 
those related to the sum of the NNLO^ term and the 
twist-six contribution [ij. Because these two latter con- 
tributions are comparable in size, they almost mutually 
cancel. For a more detailed understanding, we have com- 
piled the calculated form-factor values for the BMS DA 
at each measured point in Tabic U (last column), includ- 



ing also the aforementioned theoretical errors. Note that 
the first error originates from items (i) and (ii), whereas 
the second one stems from item (iii). 

As one clearly sees from Fig. |4j these predictions 
(shaded strip in green color), obtained with LCSRs 
within the standard collincar factorization scheme of 
QCD, cannot reproduce the significant rise of the high-Q 2 
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tail of the 7*7 — > n° BaBar data Q- But, surprisingly, 
they are in very good agreement with the BaBar data for 
the two-photon r\ and 7/ decays (HJ using the description 
of the "q — rj mixing in the quark flavor basis [42[ . These 
data are also shown in the fifth column of Table fl] Un- 
der the assumption that the state \n) = (\uu) + \dd))/\/2 
has a DA similar to that of the pion, one can link the 
7*7 — > \n) transition form factor, multiplied by 3/5, to 
the form factor 7*7 — > 7r°, where the prefactor arises 
from the quark charges. Strictly speaking, one should 
compare the state tt° not with the physical particles r\ 
and rf but with the other neutral states 771 and r/g of 
the pseudoscalar meson sector. The observed discrep- 
ancy of the two BaBar data sets indicates that the struc- 
ture of 778 is very different from that of 7r°, which looks 
implausible in view of the fact that both these states 
belong to the octet (a detailed discussion of this issue 
is given in [33| and [34|). Whether this is a true dy- 
namical effect or accidental cannot be answered rigor- 
ously at present. One realizes from this figure that two 
BaBar data points for the 7*7 — > it measurement Q 
at Q 2 = 12.7 GeV 2 and 27.3 GcV 2 also lie inside the 
shaded strip, being perfectly in line with the data for the 
\n) form factor. Moreover, our predictions (shaded strip) 
approach, together with these data, the asymptotic QCD 
limit \Z2Jt, from below. In marked contrast, the predic- 
tions derived by Agaev et al. @ — single solid lines in 
this figure — fail to comply with these data and exceed 
considerably the asymptotic QCD limit. Unfortunately, 
they even fail to provide agreement with the 7*7 — » 7r° 




1 10 100 

FIG. 4: (color online). Predictions for the scaled form factor 
Q 2 F 1 171 (Q 2 ) using a logarithmic scale for Q 2 . The predic- 
tions using the BMS "bunch" of n DAs are displayed in the 
form of a shaded (green) strip, with the solid line denoting 
the BMS model [171 ] . The width of the strip contains vari- 
ous theoretical uncertainties explained in the text. The two 
.single solid (blue) lines reproduce the predictions of Agaev 
et al. |J with their models I and III, whereas the prediction 
using their model II — not shown — would be in-between. [For 
the sake of comparison, we extended these predictions beyond 
40 GeV 2 in the form of dashed lines up to the remote point 
of 112 GeV 2 ]. The experimental data are taken from various 

experiments [IMS, El]. 












!^^ + 




— ^\ 
/ - 1 

1 V ' 

- 1 1/ 

1 --"1 

- 



0.20 0.22 0.24 0.26 0.28 0.30 0.32 

FIG. 5: (color online). Predictions for the moments (£ 2 },r 
and {£ 4 }tt at the lattice scale /j 2 at = 4 GeV 2 . The solid (blue) 
ellipse corresponds to our choice of M 2 , whereas the dashed 
(red) one results when using M 2 = 1.5 GeV 2 . The vertical 
lines show the range of values computed on the lattice: dashed 
line — [3jJ; dashed-dotted (violet) line — [iol ]. 

BaBar data as well. Indeed, as one sees from Fig. |U the 
characteristic rise of these data with Q 2 beyond 10 GeV 2 
cannot be reproduced, mainly because the enhancement 
provided by a large and positive coefficient 04 is not re- 
ally sufficient — even augmenting it with the inclusion of 
some more higher coefficients as, . . . — so that the form 
factor gradually flattens out instead of increasing with 
Q 2 . However, in view of the fluctuations of the BaBar 
data at higher Q 2 values, one cannot rule out a flatter 
behavior than the fit in Eq. (JSJ) proposed by BaBar. To 
quantify the above statements and facilitate comparisons 
with other approaches, we collect in Table [TTI the statis- 
tical properties pertaining to the data fits, using various 
sets of data and selected pion DA models. Let us also 
mention that the variation of the lower bound of Q 2 in 
the statistical analysis, e.g., the exclusion of the lowest 
six experimental points to prevent a strong influence of 
twist-4 uncertainty, does not affect the pattern of the % 2 df 
values of the BMS DA, Asy DA and the CZ DA, shown 
in Table Table [El A detailed analysis of the variation of 
the lower bound in Q 2 has been given in flU ] . 

Be that as it may, the most immediate conclusion is 
that the BaBar data for the pion and the \n) state seem 
antithetical and cannot be equally reproduced within the 
standard QCD scheme. 

The accuracy of the experimental data on 7*7 — > tt° 
and the precision of current lattice simulations in extract- 
ing constraints on the second moment (£ 2 )tt of the pion 
DA have both reached a level that can be used to deter- 
mine the range of values of the next moment where 
both experimental data and lattice computations over- 
lap. This possibility was first pointed out in (43| and the 
following range of values was extracted from the la error 
ellipse of the CLEO data Q in conjunction with the lat- 
tice constraints for (£ 2 )„ from 0: (£ 4 ) T e [0.095-^0.134] 
at M 2 at = 4 GcV 2 and for M 2 = 0.7 GcV 2 . Here 
we refine this procedure in the following way: first, 
we map the (02,04) plane onto the ((£, 2 )tt, (£ 4 )tt) plane. 
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TABLE II: Statistical properties of selected pion DA models (first column) with corresponding coefficients a n (second column) 
used in the calculation of the pion-photon transition form factor by means of LCSRs. The last two columns show the values of 
Xndi = X 2 / n df (with ndf — number of degrees of freedom) for the data in the CLEO regime and for the whole set of the data, 
respectively. 

7T DA model (q 2 , a 4 , . . x j dl (1-9 GeV 2 ) X g d f (1-40 GeV 2 ) 

BMS DA (X) (0.141,-0.089) 0.5 3.1 

Asy DA (♦) (0,0) 4.7 7.9 

CZ DA (■) (0.394,0) 32.3 25.5 

Model I [9] (0.084,0.137,0.088,0.063,0.048,0.039) 2.8 2.4 

Model III [9] (A) (0.104,0.123,0.039) 3.2 2.8 



Then, we combine the mapped la error ellipse (area en- 
closed by a solid blue line) in Fig. [5] for the data sets 
in the CLEO regime (cf. Fig. with the constraints 
from [39| and [ic| (Table XII there), the aim being to 
extract that range of values of (S, 4 )^ where they over- 
lap. The corresponding results at the typical lattice scale 
= 4 GcV 2 arCj respectively, (i) (£ 2 ) n G [0.23^-0.29] 
and (£% G [0.102 -=-0.122], (ii) G [0.26-^0.29] and 

($ 4 ) 7r G [0.11 -f- 0.122]. These windows were extracted by 
employing a Q 2 -dependent Borel parameter, as before in 
this analysis. Using instead M 2 = 1.5 GeV 2 , as in 
one would obtain only a small intersection of the cor- 
responding region, shown in Fig. [5] by the dashed (red) 
line, with the lattice constraints of [39(, amounting to the 
value (£ 4 )tt — 0.1, while there would be no intersection 
at all with the lattice estimates of [io| • The sensitivity of 
(S, 2 )^ on the choice of the Borel parameter M 2 provides 
a handle to consider the existing lattice computations 
[39l . l40j as providing independent evidence in support of 
the original prescription for the Borel parameter in the 
LCSRs for the pion-photon transition form factor [H, [ll[ . 

V. SUMMARY AND CONCLUSIONS 

In summary, three different sets of experimental data 
on the pion-photon transition form factor have been an- 
alyzed with the help of QCD LCSRs along the lines of 
our earlier works in I15l . [iff and adopting the reason- 
ing expressed in [H, lll| . The following key ingredients 
have been taken into account: (i) the NLO radiative 
correction, (ii) the twist-four contribution, and (hi) the 
NNLO^ QCD correction [H, Hi] together with the re- 
cently computed 0] twist-six contribution, these two by 
means of theoretical uncertainties. This is possible in 
our approach because in the range of the Borel param- 
eter adopted, M 2 < 1.0 GeV 2 , these two corrections, 
though with opposite signs, have almost the same abso- 
lute (not large) magnitude. Carrying out the analysis, 
evidence has been obtained that the data from CELLO 
0, CLEO @, and BaBar between 1 and 9 GeV 2 favor 
a pion distribution amplitude with endpoint suppression, 
like the BMS model [T3 (see the left panels of Fig. Hand 
Fig. [3]). The key for the overlap between the error el- 



lipses of the data in the (02,04) plane with the allowed 
region from NLC SRs is a negative value of the Gegen- 
bauer coefficient \a^\ < 02. Beyond 10 GeV 2 , the best fit 
to the BaBar data on F 1 7_>7r (Q 2 ) requires a sizeable 
coefficient a§ comparable with the lower coefficients (cf. 
right panel of Fig. [5]). 

The comparison of our theoretical predictions for the 
scaled form factor Q 2 F 1 ~ f7T (Q 2 ) in Fig.0]makcs it appar- 
ent that there is an antithetic trend between the BaBar 
data on the pion-gamma transition form factor Q and 
those extracted from the the 7*7 — > 77(77') transition form 
factors [25j]. While the latter agree very well with the 
BMS predictions (solid line in Fig. [1] and shaded strip 
in Fig. [4]), the high-Q 2 tail of the n° BaBar data re- 
quires a pion DA with sizeable (or even growing with n) 
higher Gcgenbauer coefficients a n , or alternative theo- 
retical schemes outside the standard QCD factorization 
approach, see, e.g., (2f| I27I l30l l45l l46j . Similar con- 
clusions were also drawn in [47| using Dyson-Schwingcr 
equations. 

This intriguing behavior of the BaBar data cannot be 
reconciled with QCD factorization in combination with 
LCSRs, despite the opposite claims expressed in 0. As 
one sees from Fig. 0] this is not possible for our choice 
of the Borel window M 2 « 0.7 — 0.9 GeV 2 , giving rise 
to the shaded strip, but also not for the higher values 
M 2 = 1.5 ± 0.5 GeV 2 employed in @ (single solid lines 
in the same figure). Even a large positive coefficient 04 
and a sizeable value of cannot provide sufficient re- 
inforcement of the 7*7 — > 7r° form factor to bridge the 
gap to the high-Q 2 BaBar data and reproduce their in- 
crease with Q 2 without the loss of the statistical accuracy. 
Quantitatively, this means that involving in the 3D data 
analysis the high-Q 2 tail of the BaBar data, makes the de- 
scription inevitably worse (Xndi — -0 m comparison with 
the analysis applicable to the CLEO regime (x 2 d f ~ 0.4). 
Remarkably, also the recent light-front holographic anal- 
ysis in Ref. [34j yields predictions which are incompati- 
ble with the strong rise of the BaBar data at high Q 2 , 
while being in agreement with those extracted from the 
7*7 — > r)(r)') transition form factors and rather close to 
our results in this region. 

Bottom line: New measurements, e.g., by the Belle 
Collaboration, may help resolving the controversies 
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around the BaBar data and understanding the dichotomy 
between the pion and the state |n) which indicates a 
strong violation of flavor symmetry. In addition, a precise 
lattice estimate of (£ 4 )7r could be instrumental in fixing 
the sign of the Gegenbauer coefficient 04 of the pion DA. 
Moreover, our detailed predictions in Table U in conjunc- 
tion with Fig. 2] may serve as a sort of 'reference model' 
for various purposes: theoretical and experimental. 
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